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Abstract 

A search for direct pair production of supersymmetric top squarks (ti) is presented, assuming the t\ decays into 
a top quark and the lightest supersymmetric particle, Xi, and that both top quarks decay to purely hadronic 
final states. A total of 16 (4) events are observed compared to a predicted Standard Model background of 
13.5^ 36 (4.4^ J'3) events in two signal regions based on JCdt = 4.7 fb -1 of pp collision data taken at yfs = 
7 TeV with the ATLAS detector at the LHC. An exclusion region in the i\ versus X? mass plane is evaluated: 
370 < mi < 465 GeV is excluded for to~o ~ GeV while mi = 445 GeV is excluded for m-o < 50 GeV. 
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A search for direct pair production of supersymmetric top squarks (ti) is presented, assuming 
the ti decays into a top quark and the lightest supersymmetric particle, X?, and that both top 
quarks decay to purely hadronic final states. A total of 16 (4) events are observed compared to a 
predicted Standard Model background of 13.51 3'j (4.4l ^'3) events in two signal regions based on 
fCdt = 4.7 fb -1 of pp collision data taken at = 7 TeV with the ATLAS detector at the LHC. 

An exclusion region in the ti versus X? mass plane is evaluated: 370 < mj < 465 GeV is excluded 
for m-o ~ GeV while m t -, = 445 GeV is excluded for m-o < 50 GeV. 
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The Standard Model (SM) is a successful but incom- 
plete model of particle interactions. Supersymmctry 
(SUSY) [IH3 provides an elegant solution to cancel the 
quadratic mass divergences that would accompany a SM 
Higgs boson by introducing supersymmetric partners of 
all SM particles, such as a scalar partner of the top quark 
(i). Like ti, direct tt* is produced primarily through 
gluon fusion at the Large Hadron Collider (LHC). The 
production cross section depends mostly on the mass of 
the top partner and has minimal dependence on other 
SUSY parameters (Tol - [l^ | . The LHC enables searches for 
direct it* production at higher mass scales than previous 
accelerators [l3l - l27j . The viability of SUSY as a scenario 
to stabilize the Higgs potential and to be consistent with 
electroweak naturalness 28, [29| is tested by the search 



for t below the TeV scale. 

In this Letter, we present a search for direct tt* produc- 
tion assuming ti — > tXi — > bWXi where i\ is the lightest t 
eigenstate and Xi represents the lightest supersymmetric 
particle (LSP) in i?-parity conserving models Hail. We 
consider events where both W bosons decay hadronically, 
yielding a final state with six high transverse momentum 
(pt) jets from the all-hadronic ti final state and large 
missing transverse momentum (E™ lss ) from the LSPs. 
The kinematics of both top quarks are therefore fully 
specified by the visible decay products. Additionally, SM 
backgrounds from all-hadronic ti are suppressed as there 
is no intrinsic _E™ 1SS except from semi-leptonic c- and b- 
quark decays. The dominant background consists of ti 
events that contain a W — > iv decay where the lepton 
(£), often of r flavor, is either lost or mis-identified as a 
jet and which have large E™ lss from the neutrino. 

The data sample was acquired during 2011 in LHC pp 
collisions at a center-of-mass energy of 7 TeV with the 
ATLAS detector 35[ , which consists of tracking detectors 
surrounded by a 2 T superconducting solenoid, calorime- 
ters, and a muon spectrometer in a toroidal magnetic 
field. The high-granularity calorimeter system, with ac- 
ceptance covering \rj\ < 4.9 [36(, is composed of liquid 



argon with lead, copper, or tungsten absorbers and scin- 



tillator tiles with steel absorbers. This data set, com- 
posed of events with a high-px jet and large E™ lss as 
selected by the trigger system, corresponds to an inte- 
grated luminosity of 4.7 fb _1 with a relative uncertainty 
of 3.9% ©ll^. 

Jets are constructed from three-dimensional clusters 
of calorimeter cells using the anti-fct algorithm with a 
distance parameter of 0.4 [39. 40]. Jet energies are cor- 
rected [4l| for losses in material in front of the active 
calorimeter layers, detector inhomogeneities, the non- 
compensating nature of the calorimeter, and the impact 
of multiple overlapping pp interactions. These correc- 
tions are derived from test beam, cosmic ray , and pp 
collision data, and from a detailed Geant4 [42( detector 
simulation [43j . Jets containing a b- hadron are identified 
with an algorithm exploiting both the impact parame- 
ter and secondary vertex information [44, 45 1. A factor 
correcting for the slight differences in the 6-tagging effi- 
ciency between data and the Geant4 simulation is applied 
to each jet in the simulation. The 6-jets are restricted 
to the fiducial region of the tracker, \t]\ < 2.5. Non-ii 
backgrounds are minimized by requiring either > 1 6-jets 
with a selection corresponding to a 60% efficiency with 
a low < 0.2% mis-identification rate (tight), or > 2 b- 
jets each with 75% efficiency but a higher ss 1.7% mis- 
identification rate per 6-jet (loose). 

The E™ lss is the magnitude of p™ ISS , the negative vec- 
tor sum of the pt of the clusters of calorimeter cells, cal- 
ibrated according to their associated reconstructed ob- 
ject (e.g., jets and electrons), and of the pt of muons 
above 10 GeV within |r;| < 2.4. Events containing 
fake E™ lss induced by jets associated with calorimeter 
noise or non-collision backgrounds (4fj|, or by cosmic- 
ray muons or poorly reconstructed muons [471 l48j |. are 
removed from consideration. Large p™ ISS colinear with 
a high-px jet could indicate a significant fluctuation in 
the reconstructed jet energy or the presence of a semi- 
lcptonic c- or 6-quark decay. Therefore, the difference in 
azimuthal angle (A^>) between the p™ ISS and any of the 
three highest-px jets in the event, A0(p T llss , jet), is re- 
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quired to be > 7r/5 radians. Fake E™ laa is also suppressed 
by requiring that the A<p between the above computed 
p™ ss and one calculated with the tracking system, using 
tracks having px > 0.5 GeV, is < 7r/3 radians. 

Events are required to have at least one jet with 
p T > 130 GeV in \rj\ < 2.8 and E™ iss > 150 GeV to 
ensure full efficiency of the trigger. At least five other 
jets having px > 30 GeV and \rj\ < 2.8 must be present. 
In addition to the jet and E™* ss requirements, events 
containing "loose" electrons [HI with pt > 20 GeV 
and \r)\ < 2.47 that do not overlap with any jet within 
AR < 0.4, where AR = ^/(Arj) 2 + (A(t>) 2 , are rejected. 
Similarly, events with muons [47|,[51| having px > 10 GeV 
and \r]\ < 2.4 that are separated by AR > 0.4 from 
the nearest jet are rejected. A jet with 1-4 tracks and 
A0(p T 1ISS , jet) < 7r/5 indicates a likely W — > rv decay. 
Events with such r-likc jets that have transverse mass 
m T = ^2p T E™ iss {l-cosA(f>) < 100 GeV are rejected. 

The presence of high-px top quarks that decay through 
t —> bW — ^ bjj in the tit* final state is exploited 
to further reduce SM backgrounds by only considering 
events with reconstructed three-jet invariant masses con- 
sistent with the top-quark mass (m t ). A clustering tech- 
nique resolves the combinatorics associated with high- 
multiplicity jet events. The three closest jets in the 
rj — <p plane are combined to form one triplet; a sec- 
ond triplet is formed from the remaining jets by repeat- 
ing the procedure. The resulting three-jet mass (rrijjj) 
spectrum is shown in Fig. [T] for the control region con- 
structed from £+jets events (defined below). There is a 
clear peak associated with the hadronically-decaying top 
quarks above a small non-ti background. A requirement 



of 80 < 



< 270 GeV is placed on each reconstructed 



triplet in the event. The kinematics of the t — > bW — > blv 
decay is also exploited to further reduce the dominant 
£+jets ti background, as the mx distribution of the p T 1ISS 
and &-jet (m^ ) has an endpoint at m t . When there are 
> 2 loose 6-jets, the for the 6-jet closest to the p T nss 
is required to be > 175 GeV. The largest m^*, calcu- 
lated for each of the four highest-px jets, is required to 
be > 175 GeV in the case of only one tight 6-jet. 

Two signal regions (SR) arc defined including the 
above kinematic and mass requirements. The first, which 
requires E™ lss > 150 GeV (SRA), is optimized for low 
toj , while the second, requiring E™ lss > 260 GeV 
(SRB), is used for higher mj . Using these signal re- 
gions, the search is most sensitive to t\t* production 
with 350 < mt < 500 GeV and m-o < mi. Sig- 
nal events are simulated using Herwig+- 1- [52| with the 
MRST2007LO* parton distribution functions (PDF) 
generated with the i\ and Xi masses at fixed values in 
a grid with 50 GeV spacing. The mixing parameters 
are chosen such that the lightest scalar top is mostly 
the partner of the right-handed top quark. The branch- 
ing fraction of ti —> tXi is set to 100% 54 1. Signal 
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FIG. 1. Three-jet invariant mass distribution of the hadronic 
top-quark candidate in the control region constructed from 
^+jets events. Data are indicated by points; shaded his- 
tograms represent contributions from several SM sources 
(with ti scaled by 0.66). The hatched error bars indicate 
the total statistical uncertainty on the expected background. 
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expectation in this control region is overlaid. 



cross sections are calculated to next-to-leading order in 
the strong coupling constant, including the resummation 
of soft gluon emission at next-to-leading-logarithmic ac- 
curacy (NLO+NLL) [IM3- 

The nominal production 
cross section and associated uncertainty are taken from 
an envelope of cross section predictions using different 
PDF sets and factorization and renormalization scales, 
as described in Ref. [Hij . The tit* cross section for 
m h = 400 GeV is rj flt -» = 0.21 ± 0.03 pb. 

In the signal region, the dominant source of SM back- 
ground is tt — > r+jets events where the r lepton is re- 
constructed as a jet. Additional, smaller, backgrounds 
include other ti — > £+jets final states, ti + V where V 
represents a W or Z boson, single top quark produc- 
tion, V+jets, and W+jets. The ti events are produced 
with ALPGEN [HI] usingthe CTEQ6L1 PDF [H^ and 
interfaced to HERWIG [H, for particle production 
and JIMMY (6(| for the underlying event model. Addi- 
tional ti samples generated with MC@NLO [i3| and 
AcerMC [H], interfaced to HERWIG and JIMMY, are 
used to estimate event generator systematic uncertain- 
ties. Samples of ti+ V are produced with MadGraph 64 1 



interfaced to PYTHIA 
generated with MC@NLO 



65 , 66J . Single top events are 
ialHI and AcerMC. The as- 
sociated production of W and Z bosons and light and 
heavy-flavor jets is simulated using ALPGEN; diboson 
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FIG. 2. The rnx distribution for r-like jets with the selection 
described in the text. Data are indicated by points; shaded 
histograms represent contributions from several SM sources 
(with tt scaled by 0.66). The hatched error bar indicates 
the systematic uncertainty on the total expected background. 
The expected signal distribution for = 400 GeV,m^o = 
1 GeV is overlaid. 



FIG. 3. The _B™ 1SS distribution in data compared to the 
SM expectation for signal region A. The hatched error bar 
indicates the systematic uncertainty on the total expected 
background. The expected signal distribution for m%. = 
400 GeV, m-o = 1 GeV is overlaid. The SM background dis- 
tributions do not include the 0.2 ± 0.2 events of all-hadronic 
tt and QCD multi-jets estimated from data. 



production is simulated with SHERPA [69] . 

All samples are passed through the Geant4 simula- 
tion of the ATLAS detector, and are reconstructed in 
the same manner as the data. The simulation includes 
the effect of multiple pp interactions and is weighted to 
reproduce the observed distribution of the number of in- 
teractions per bunch crossing. SM event samples are nor- 
malized to the results of higher-order calculations using 
the cross sections cited in Ref. [7(| except for the £+jets tt 
background. This sample is renormalized by a factor that 
scales the tt expectation to agree with the observed data 
in a control region (CR) kinematically close to the signal 
region but with little expected signal. The CR is con- 
structed from events containing one muon or one "tight" 
electron [49j with p^ > 30 GeV consistent with originat- 
ing from a PF-boson decay (40 < m^. < 120 GeV) and > 5 
jets, where is the transverse mass of the electron or 
muon and p!j? lss . The lepton must be isolated such that 
the scalar sum of tracks in a cone of AR < 0.2 around 
the lepton, excluding the track of the lepton, is < 1.8 GeV 
for the muon or is < 10% of the electron pr, respectively. 
The jet, 6-jet, and E™ lss requirements remain the same as 
the standard signal selection; however, some topological 
constraints are relaxed (A0(p T llss , jet) > tt/10 radians 
and rrijjj < 600 GeV) and others removed (m^ ) to gain 
statistics. The tt purity in the control region is > 80%; 
the expected signal contamination is < 3%. The lepton is 



treated as a jet of the same energy and momentum, mim- 
icking the effect of the r lepton. Effects of the additional 
fimiss f rom ^ ue T neu t r ino are smaller than the statisti- 
cal uncertainties. The scale factor needed to bring the 
> 6 jet ^+jcts ALPGEN tt events into agreement with 
the data after recalculating all quantities except i?™ lss 
is 0.66 ± 0.05; the uncertainty quoted here is statistical 
only. This scale factor is used in Figs. UHll The nor- 
malization is validated with an orthogonal if-dominated 
sample created from SRA by selecting events with r-like 
jets; the requirement on is removed to increase the 
sample size. The tot of r-like jets is shown in Fig. [21 
where the tt sample has been normalized as described 
above. Expectations from the simulation agree with the 
data within uncertainties. Contributions to the signal re- 
gion from QCD multi-jet and all-hadronic tt production 
are estimated with a data-driven technique 7l| . Jets are 



smeared in a low-E^j? lss data sample using response func- 
tions derived from control regions dominated by multi-jet 
events. Only 0.2 ± 0.2 such events remain in SRA after 
the full event selection. 

The E^ iss distribution in SRA is shown in Fig. [3] for 
data, for the SM backgrounds, and for expectations of 
t\t\ production with = 400 and m-o = 1 GeV. Num- 
bers of events and combined statistical and systematic 
uncertainties, for both SRA and SRB, are tabulated in 
Table[U Uncertainties in the event generators, including 
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TABLE I. The numbers of expected events for the SM back- 
grounds and an example SUSY signal point, and the observed 
number of events in data. The 95% CL S upper limit on the 
observed (expected) visible cross section is appended below. 



SRA 
> 150 GeV 



SRB 
> 260 GeV 



IT 


9.2 


±2.7 


2.3 


± 


0.6 


ti + W/Z 


0.8 


±0.2 


0.4 


± 


0.1 


Single top 


0.7 


±0.4 


0.2 


+ 


0.3 
0.2 


Z+jets 


1.3 


+ i.i 

- 1.0 


0.9 


+ 


0.8 
0.7 


W+jets 


1.2 


+ 1.4 
- 1.0 


0.5 


± 


0.4 


Diboson 


0.1 


+ 0.2 

- 0.1 


0.1 


+ 


0.2 
0.1 


Multi-jets 


0.2 


±0.2 


0.01 


± 


0.02 


Total SM 


13.5 


+ 3.7 
- 3.6 


4.4 


+ 


1.7 

1.3 


SUSY (m^m-o) = (400,1) 
Data (observed) 


GeV 14.8 


±4.0 


8.9 


± 


3.1 




16 
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the impact of initial- and final-state radiation, are the 
dominant source of systematic uncertainty of 28% (23%) 
in SRA (SRB). Other major sources of uncertainty in- 
clude 22% (32%) for the jet energy calibration, 6.5% 
(6.8%) for jet energy resolution, 5.9% (6.2%) for 6-jct 
identification, and 1.4% (1.5%) for E™ iss in SRA (SRB). 

The number of observed events in the data is well 
matched by the SM background. These results are in- 
terpreted as exclusion limits for and m-o using a 
CL S likelihood ratio combining Poisson probabilities for 
signal and background [72[ . Systematic uncertainties are 
treated as nuisance parameters assuming Gaussian distri- 
butions. Uncertainties associated with jets, &-jets, E!f lss , 
and luminosity are fully correlated between signal and 
background; the others are assumed to be uncorrelated. 
The expected limits for the signal regions are evaluated 
for each (m^, m~o) point; the SR with the better ex- 
pected sensitivity is used for that point. The expected 
and observed 95% C.L. exclusion limits are displayed in 
Fig. |4] Top squark masses between 370 and 465 GeV are 
excluded for m~o ~ GeV while m; = 445 GeV is ex- 

Xi L l 

eluded for m-o < 50 GeV. These values are derived from 

i 

the —la observed limit contour to account for theoreti- 
cal uncertainties on the SUSY cross sections. The 95% 
CL S upper limit on the number of events beyond the SM 
in each signal region, divided by the integrated luminos- 
ity, yields limits on the observed (expected) visible cross 
sections of 2.9 (2.5) fb in SRA and 1.3 (1.3) fb in SRB. 

In conclusion, we have presented a search for the direct 
production of tit\ in the all-hadronic tXitXi decay chan- 
nel, assuming B{h -» tXi) = 100%. A total of 16 (4) 
events are observed compared to a predicted Standard 
Model background of 13. 5± f 7 6 (4.4! events in two 
signal regions based on JCdt — 4.7 fb -1 of pp collision 



FIG. 4. Expected and observed 95% CL S exclusion limits in 
the plane of m-o vs. m,f , assuming B(ti — > tXl) = 100%. 
The dashed line shows the expected limit at 95% C.L. with 
the shaded region indicating the ±1<t exclusions due to exper- 
imental uncertainties. Observed limits are indicated by the 
solid contour (nominal) and the dotted contours (obtained by 
varying the SUSY cross section by the theoretical uncertain- 
ties). The inner dotted contour indicates the excluded region. 



data taken at y/s = 7 TeV with the ATLAS detector at 
the LHC. No evidence for tit* is observed in data and 
95% C.L. limits are set on t\t\ production as a function 
of mz and m^o . 
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